Quantifying mixed-state quantum entanglement by optimal entanglement witness 
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We develop an approach of quantifying entanglement in mixed quantum states by the optimal 
entanglement witness operator. We identify the convex set of mixed states for which a single witness 
provides the exact value of an entanglement measure, and show that the convexity and symmetries 
of entanglement or of a target state considerably fix the form of the optimal witness. This greatly 
reduces difficulty in computing and experimentally determining entanglement measures. As an 
example, we show how to experimentally quantify bound entanglement in four-qubit noisy Smolin 
states and three-qubit Greenberger-Horne-Zeilinger (GHZ) entanglement under white noise. 

PACS numbers: 03.67.Mn, 03.65.Ud 



Introduction. ~ Entanglement is at the center of foun- 
dational issues about quantum strangeness such as non- 
locality. It is a key resource of quantum technologies [1] . 

A long-standing issue is how to detect and quantify 
entanglement [T] [2] . Entanglement measures have been 
proposed, to quantify entanglement, to estimate its role 
in quantum information tasks [TJ [5], and to character- 
ize its various aspects of such as multipartite entangle- 
ment classes [2] and correlations in many-body states and 
phase transitions [51. Many useful measures £, such as 
entanglement of formation, concurrence, 3-tangle, etc., 
are defined to be computable for pure states ^, and ex- 
tended to mixed states p via the convex roof construction 



f(p)= inf Vp,£(|V.)), 
{pi,il'i} '—' 



(1) 



where the infimum is taken over all possible pure-state 
decompositions p = Yl,iVMt)ki'i\ with YjiVi = 1, Pj > 0, 
and (V'ilV'i) = 1- This construction satisfies the require- 
ment T] that a convex measure E does not increase under 
local operations and classical communications. 

It is hard to compute and to experimentally determine 
E for mixed states. For two qubits, concurrence is com- 
putable [3] and experimentally determinable [51 [3] . For 
higher dimensional or multipartite cases, however, the 
computation of f (p) requires an impractical numerical 
task of exploring all possible pure-state decompositions 
of p, except for rare specific states ffHIO] with very low 
rank or high symmetry. This difficulty obstructs quanti- 
tative researches on the features of mixed-state entangle- 
ment such as nonlocality, fragility under noise, dynamics, 
and applicability to quantum technologies. 

In this work, we develop a quantification approach, 
based on a physical observable known as the optimal en- 
tanglement witness operator [TTJ [T^] whose expectation 
value provides the exact value of an entanglement mea- 
sure for a state; the optimality is differently defined in 
Ref. [13]. We reveal the basic relation (see Theorems) 
between the optimal witness Xp for f (p) and the opti- 
mal pure-state decomposition of p in Eq. ([l]). It allows 
one to identify the convex set of mixed states for which 



a single witness provides the exact value of a measure, 
and combining with symmetries of entanglement (e.g., 
SLOCC invariance) and a target state, it considerably 
fixes the form of the optimal witness; SLOCC stands for 
stochastic local operations and classical communications. 
These findings greatly reduce the difficulty of computing 
and experimentally determining E, and apply to general 
convex measures. As an example, we show how to exper- 
imentally quantify bound entanglement in 4-qubit noisy 
Smolin state [131 [12] and 3-qubit GHZ entanglement [TB] 
under noise such as full-rank white noise. 

Approach. - Our approach for entanglement quantifica- 
tion starts with an equivalent expression oi E{p) [111117) . 



E{p)= sup Tv{Xp), 

xeM-H 



(2) 



I« = {X\TT{XTr^) < Eim for ^|V') G H}. 



TT^ = |V')(V'| is the projector on a pure state ^, Tr(-) is the 
trace, and M^ includes the null operator to have E > 0. 
A physical operator X G M.-^ does not overestimate 
E{\tp)) for all pure states V' in a Hilbert space H, hence 
—X is a witness [12 \T^ [TB] detecting entanglement; see 
Fig. [IJa). Among X's, the optimal one with the largest 
expectation value of Tr(Xp) detects E{p), and a less opti- 
mal one provides a lower bound oiE{p) [TTIIT^ITTIITM^ . 
The optimal witness exists as a support hyperplane of a 
convex set of states [see Fig. [l](b)], and is an observable 
in principle experimentally accessible [^. Hereafter, we 
would call X a witness, although —X is a conventional 
one, and accordingly redefine the optimal one: 

Definition. The optimal witness Xp G M-^ is defined 
relative to p and E. It satisfies Ti{Xpp) — E{p). 

To compute E{p), one needs to find Xp. We below show 
the basic properties of Xp useful for finding Xp efficiently. 

Theorem I. The optimal witness Xp for p is also opti- 
mal for all the pure states ip°^ of the optimal decompo- 
sition of p giving E{p) = J2iP7^i\'^7))- 
Proof by contradiction. Suppose that Xp is not the op- 
timal witness for ^ip°^, Tr(Xp7r^°p) < E{^j^). Then, 
TriXpp) = E.PrTr(X,7r^op) < t.^ff (I^D) = ^(p), 
contradicting the fact that Xp is optimal for p. D 



Tr{Xp) = 



Tr(Xpp) = e(p) 




Tr(Xp,) = <^ S{p2) 
U(P3) 




FIG. 1: (Color Online) Optimal entanglement witness, (a) A witness operator X detects entanglement in a quantum state p, 
as it divides the set (shade area) of states into two parts with and without separable states by the sign of the expectation value, 
(b) Hierarchy of convex sets S{x), each consisting of the states p' with an entanglement measure £{p') < x; darker for smaller x. 
The witness Xp optimal for p quantifies entanglement by £{p) = Tr(Xpp), which corresponds to the support hyperplane (solid 
line) of S(a; = £^(p)) at p. (c) When X is optimal for {p^} (i = 1, 2, • • • ), so is it for all states in the convex hull, conv({pi}), of 
{pi} (shade triangle). For ^ p € conv({pi}), one computes E(j)) by Tr(Xp), without exploring pure-state decompositions of p. 



Theorem 2. Consider the set ¥x C "H of the pure 
states llV'i)} for which a witness X E M-u is optimal, 
Tt{Xtt^.) — SQipi)). Then, for any convex mixture 
p — J2iPi''^ipi G conv(Px), X is optimal, Xp = X, and 
'^iPi''^ipi is the optimal decomposition of p. 
Proof. £{p) > Tt{Xp) = E^P^MX^^i,,) from Eq. (2), 
^(P) < E^P^£i\^^)) from Eq. (1), and Tr(X7r^J = 
£{\^,)). .-. Sip) = E.P.^dV-.)) = Tr(Xp). D 

The theorems connect Xp and the optimal pure-state 
decomposition of p in Eq. (fTl) . They are valid for all con- 
vex measures. Crucial are their consequences [Fig. IIJc)]: 
One can check the optimality of a witness X for p, by see- 
ing whether p £ conv(Pjc). Namely, one computes £{p) 
or its lower bound, by optimizing the form of Xp, or by 
guessing the form and checking its optimality for p (with 
possibly avoiding heavy computation). Moreover, with a 
single Xp, one can obtain the analytic expression of £ or 
experimentally determine £ for all p G conv(Px ). This 
considerably reduces not only computational but also ex- 
perimental efforts. For example, when a target state p 
is affected by noise in experiments, the optimal witness 
Xp of p can still give the exact value of an entanglement 
measure of the affected state or a faithful lower bound, if 
the affected state is belonging to or close to conv(Pxp). 

We further provide the restrictions on the form of Xp 
by the range Tl{p) and rank Tp of p. The Hilbert space H 
for M^ is TZ{p) or the full Hilbert space Hf of the system. 
Theorem 1 ensures another restriction useful for large rp-. 

Corollary. The number of linearly independent pure 
states in Vx should be larger than or equals to rp. 

Symmetries of £ or of p also restrict Xp. Many useful 
measures characterize SLOCC invariant entanglement. 
In this case, all pure states with finite £ are connected, 
by SLOCC operations (tensor products of local opera- 
tors with determinant 1), to a maximally entangled state. 
This connection simplifies Xp-, see, e.g., Eq. ^. 

Moreover, when p has (higher) symmetries, it is enough 
to consider the (simpler) form of Xp symmetrized by 
the same symmetries, as Tr(Xpp) indicates. For highly 



symmetric states, our approach reproduces (hence cov- 
ers) previous theoretical results [THTO]. and also enables 
experimental quantification, contrary to the previous 
works; for example, for an isotropic state and a bipar- 
tite Werner state, Xp has only one parameter, hence one 
analytically obtains Xp and £ easily. For states having 
low symmetry, multipartite, high rank, and/or nontrivial 
decomposition of Eq. (fTj) , our approach is more useful for 
obtaining Xp and £ than the previous works; see below. 
Noisy Smolin states. -For illustration, we first quantify 
entanglement in four-qubit noisy Smolin states 



Pm{p) = (l-p)ps- 
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by geometric measure of entanglement j9j £g; ctj's are 
Pauli matrices and / is the identity. For p £ [0,2/3), 
P^sip) h&s bound entanglement experimentally more re- 
liable than Smolin state ps EH US] ■ It was realized ex- 
perimentally |15| . Its entanglement has never been quan- 
tified, while £q{ps) was computed in Ref [H]. 

We derive a general condition of X for £q, which 
greatly reduces computation costs. With the set U of 
separable pure states, £g is defined [S] as £'g(|V')) = 
1 — max|s)gi[j |(V'|s)|^ for pure states ip, and extended to 
mixed cases via Eq. ([l]). Since (f/'l^ [■(/') < £g(|V')) < 
l-|(V'|s)|^ {^p\X + Tr,\iP) < Ifor^ls) G U and ^X. The 
equality holds, when X is optimal for ?/) and \s) is the 
state in U with maximal |(V'|s)|^. Since any X is optimal 
for at least one pure state (see Corollary), it satisfies 



maxAmax[^ + 7rs] = 1, 

\s)GV 



(3) 



where Aniax[^] denotes the maximum eigenvalue of A. 



We now turn to pNs(p)- 



Its symmetries restrict 
Xp^^ [p] to the form of X{p) ^ 4a(p)ps + P{p){I - 4ps) 
with only two real parameters a > and /3 < 0. From 
the witness condition of Eq. (Is]), one has a~^ -\- P~^ = 2 



or a ~ /3 = 0. By maximizing Tr{X p^s{p)) = [(4 - 
3p)a + 3p/3]/A over a and /3, we obtain for p G [0, 2/3) 



Scipmip)) 



v/3p(4-3p)]/4, 



(4) 



"(P) = [1 - yS^l/a, and /?(p) = [1 - ^/^]/2, while 
fG(pNs(p)) = ab) = Pip) = ioT p e [2/3,1]. The 
analytic result of £'g(pns) agrees with previous findings 
of separability for p e [2/3, 1] [14 and Soips) = 1/2 0, 
and does not require heavy optimization, which might be 
necessary in previous approaches [TJ [H] • 

Xpj^g is accessible in experiments. The construction 
of the common component of Xp^g, X]i ^f ' which was 
done in Ref. [TF, provides all Xpj^g(p)'s. 

3-qubit GHZ entanglement.- We next consider three- 
qubit GHZ entanglement. Using its SLOCC invariance, 
we derive the general form of the optimal witness. 

In three qubits, there are two types of genuine tri- 
partite entanglement, GHZ and W PM2S]- Their rep- 
resentative states are |GHZ) = 4|(|000) + (111)) and 
|W) = ^(jlOO) + |010) + |001)). To quantify GHZ en- 
tanglement, we define a new entanglement measure Ts [6] 
by TadV')) = v^Tad^)) for pure states, and extend it to 
mixed states via Eq. (fTl) , where t^ is three-tangle [17] . It 
is "extensive", Tad^)) = (V'jV') %{\ip) / ^/W^)) . result- 
ing in the property useful for Xp that Ts is SLOCC invari- 
ant for pure and mixed states, T^ip) — Tz{OpO^). Any 
pure state in GHZ\W class is transformed into |GHZ), 



1^) = y^sM) Ov,|GHZ) = ^/mW) O,0|GHZ), (5) 

by an SLOCC operator O^ [21]. We call Ta extensive 
three-tangle, and use it instead of three-tangle T3, since 
T3 is not SLOCC invariant for mixed states [50] . 

We demonstrate how to construct Xp for Ta. Theorem 
1 ensures that Xp should be optimal also for at least one 
pure state V in GHZ\W class, Tr(Xp|i/;)(7/;|) = TadV-))- 
This property and Eq. ^ resuU in Oj,XpO^ |GHZ) = 
I GHZ) , which gives the general form of Xp [30] , 

Xp = (TTGHZ+H-zi/) 0V(1-A*), O e SLOCC, (6) 

where ttghz = | GHZ) (GHZ |. It generalizes the widely- 
used witness of ttghz — 3//4 |2B|, and indicates that 
for p with finite Ta, p is optimally decomposed into 
one GHZ\W-class state and other W-class states, /i — 
maXi^gw('/'|7i'GHZ+H|(/3) prevents Xp from overestimating 
Ta for all pure states in %. To obtain p, it is enough to 
consider pure states ^p (with (i^|(p) = 1) in W class [50] . 
H is a sum of pure-state projectors satisfying < H < 
pi and Tr(H7rGHz) — 0. For states with rp > 2, only 
specific forms of H satisfy the Corollary. For example, 
H = Attw = A|W)(W| cannot be used for rp > 4, since 
¥x contains only four states. Our numerical search [30] 
might imply that the Corollary is satisfied only if H is 
invariant under some of the symmetries of ttghz and / 




FIG. 2: (Color Online) Three-qubit mixed-state GHZ entan- 
glement. The exact value of GHZ entanglement measure Ts 
is computed for the 3-qubit full-rank state PGGi(P''i')' ^y us- 
ing the optimal witness operator Xp for Pqqj(p, g). The 
physical region of {p, q) is defined by < p, g, (1 — p — g) < 1. 



such as permutation, exchange of qubit index, 0-1 flip, 
local phase rotation of C/(a, /3) = ( q g?a ) <8 ( q g?;? ) <8 
(J ^_i(^+^) ) with real a and /3, etc. It will be valuable 
to prove this conjecture. A more symmetric form of H 
has bigger P^ and smaller number of parameters. For 
instance, Hg = >^m¥L^GWL+T.^jkf 000 ,111 ^ijfcT^^jfe has the 
biggest set of Px , as it is invariant under all available 
local phase rotations, where [GHZ) = (|000) - |lll))/%/2 
and TTijk — \ijk){ijk\. Symmetric forms such as Hg are 
useful for the full-rank cases of rp = 8. 

We discuss the efhciency of our strategy for Ta. For a 
simple form of p, one guesses a trial witness form X, op- 
timizes it, checks the optimality whether p G conv(Px), 
then computes £{p) — Tr(A'p). This procedure is use- 
ful for p with symmetries, for which one chooses X with 
the same symmetries. On the other hand, for a complex 
form of p, one fully optimizes the form of Xp in Eq. ([6| . 
This optimization still has a cost much cheaper than the 
direct pure-state decomposition of Eq. ([l]). For rp — 8, 
Xp has 72 optimization parameters (48 for deciding the 
eigenstatcs of H, 6 for the eigenvalues of H, and 18 for O), 
while the direct decomposition has hundreds to a thou- 
sand (roughly 2r^) of parameters ;29,. If the conjecture 
about the symmetric forms of H is true, the parameter 
number of Xp is further reduced to at most 40; Xp with 
Hg has 24 parameters. Even in the case that the con- 
jecture is false, a symmetric H is useful, as it gives large 
Pjf and at least a good lower bound of Ta . 

We provide examples. We construct Xp for arbitrary 
mixtures of |GHZ), |GHZ) and the white noise 7/8, 
Pggi(P'9) = (1 - P - 9)7i'GHZ + ptTqu^ + 9T/8, 



"'^Pggi "" 



(ttghz + Att^jj2 - pl)/{l - p)- 



(7) 



Here, the symmetries of /Cqqj fix H = AtTqjj^- Checking 
the optimality, we fix A, and compute Ta in Fig. [2] 



For any mixture pQi{q)[= Pqqj(p = 0, q)] of |GHZ) and 
1/8 with (7 < go - 0.304, Xp^, = Xp^^^{fi ~ 0.750, A ~ 
0.433) is optimal, independent of q. The result of Ta is 

7i = max[0, Tr{Xp^^pGi)] = 1 - q/qo for g < qo, (8) 

Ts = for g > go — 0.304. The optimal decomposi- 
tion is PGiiq) = 79(g)7rGHZ + (1 - %{q))pz for q < go, 
and pGi(g) = T^Pz + f^l;! for g > go, where pz = 
(1 - go)7rGHZ + go^/8 oc J^i'^Zi is decomposed by W- 
class states |Zi) with Tr(XpQj7rzJ = 0; |Zi)'s are given 
in Ref. |30]. This computation has a cost significantly 
cheaper than the direct decomposition of Eq. (II]) which 
has 239 optimization parameters for pci- 

We consider another state, the rank-2 mixture 
PGw{p) = (1 -p)7i'GHZ +P7rw, < p < 1. Xp depends on 
the choice of "H between TZ{pGw) and Hi (the full space). 
For H = 7?.(pgw), a single witness quantifies Tsipcw), 
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xilL = 



PGW 



TI'GHZ 



(Pq ^ - l)'^w, 



Ts - msi40,TriX('>pGw)] 



PGw' 



(9) 
1 -p/Po ior p<po, 



Ts — for p > po — 0.373; the optimal decompo- 
sition of pgw(p) is given in Ref. [30]. For H = Hf, 

XpGw = liniA^i ihj.i'^GHZ + Attw - pi) gives 7^(pgw) 
asymptotically with lim>_yi. 

The examples are useful for understanding the fragility 
of GHZ entanglement under noise. Pggi includes the 
white noise and the dephasing of the relative phase 
between |000) and |111). GHZ entanglement in p = 
(1 — a)7rGHZ + oH' decreases as a increases, and van- 
ishes at a = 0.5 for H' = tt^jj^, 0.373 for ttw, and 0.304 
for 1/8; for tt ^^^^ , GHZ entanglement revives and in- 
creases with a > 0.5. The white noise destroys GHZ en- 
tanglement more than tt q^^ and ttw, and more severely 
for states with more qubits; a — 0.304 is smaller than 
/3 = 2/3 at which Bell-state (|Bell)) entanglement in 
P2 = (1 - /3)|Bell)(BeU| + I3I/4 vanishes. 

From our finding, one can quantify GHZ entanglement 
in experiments. When |GHZ) is prepared, it normally 
becomes Pqgi "^^^ ^^ noise. Assuming that the prepared 
state pcxp has the form of Pggi (p^ i) ' '^^^ estimates Ts by 
Xp : One measures ttghz and tTqjj^ Pi, and computes 
the largest value of Tr(Xp^^jPoxp) with varying (/i. A) 
over the space where Xp is a witness. Then, it is the 
exact value (a faithful lower bound) of Ta if the assump- 
tion is correct (incorrect). This procedure is powerful, as 
Ta is obtained from minimal information about poxp- 

Conclusion.- Our approach of optimal witness has 
great advantage over previous methods of optimizing 
state decomposition. It offers a simple way of theoreti- 
cal and experimental quantification of entanglement pre- 
pared in laboratories (which is usually in a simple state 
such as pns and Pggi)' ^^^ stimulates researches on mul- 
tipartite or high-dimensional mixed-state entanglement. 
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